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Abstract 

We develop a theory of the EPR-like effects due to neutrino oscillations in 
the 7r — > \xv decays. Its experimental implications are space-time correlations 
of neutrino and muon when they are both detected, while the pion decay 
point is not fixed. However, the more radical possibility of /i-oscillations in 



experiments where only muons are detected (as it has been suggested in [hep- 
|ph/9509261| ), is ruled out. We start by discussing decays of monochromatic 
pions, and point out a few "paradoxes" . Then we consider pion wave packets, 
solve the "paradoxes", and show that the formulas for \iv correlations can 
be transformed into the usual expressions, describing neutrino oscillations, 
as soon as the pion decay point is fixed. 



1 Introduction 

It was suggested recently [I], |2| that the much searched for hypothetical neu- 
trino flavour oscillations can cause space-time oscillations of the observation 
rate of the beam of muons from the 7r-meson decay. As emphasized in [IB], 
the same claim applies to other charged leptons, in particular to r-leptons 
in the decay W — ► tv t . If true, this phenomenon could provide a powerful 
experimental method for the (indirect) observation of neutrino oscillations. 
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We think that argumentation of [|1|] is unfortunately erroneous, and there 
are no oscillations of charged leptons if neutrinos are not observed. However, 
in coincidence experiments, in which both the charged lepton and neutrino, 
born in the same decay, are detected, specific EPR-like oscillating correla- 
tions || can show up (see Section 3.2). Once neutrino oscillations [[f]-|§ 
per se are discovered, this could be of interest for the next generation of 
experiments. 

If the neutrino mass matrix is not flavour diagonal then the "current" 
neutrinos v a (a = e, /i, r) are non-trivial mixtures of the mass eigenstates, v n 
(n = 1, 2, 3) with masses m\ > rri2 > m^. 

^a = J2 U ^n, (1) 

n 

i.e. the relevant part of the Standard Model Lagrangian is 

C u = gJ2 hWu a + H.c = gJ2 U an l a Wu n + H.c. (2) 

a a,n 

In what follows we consider a "toy model" with only two charged leptons, 
muon and electron, and two neutrinos, V\ and v^- 

L u = g (pW{pi cos 9 + u 2 sin 9) + eW{-v x sin 9 + v 2 cos #)) + H.c. (3) 

Also, we neglect the width of the emerging muon (as the oscillations discussed 
in ref. JU survive in the limit of stable /z; in any case all the widths are easy 
to restore in the formulas). 

We consider the decay ti — >• /zz/, and analyse three types of experiments: 

A) when both neutrino and muon are detected; this is an experiment with 
two detectors in coincidence; 

B) when only muons are detected; 

C) when only neutrinos are detected. 

Evidently, B and C are single-detector experiments, and B is considerably 
easier than C and A. 

Our analysis will demonstrate that the probability to detect both fi and 
v (case A) oscillates as a function of the distance d between the muon and 
neutrino detection points and of the time interval At between "clicks" of 
the two detectors. Moreover, oscillation length and oscillation frequency 
are different from the standard values L = 2E v /(m\ — ml) and L^ 1 , 
respectively. 
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These oscillations disappear completely in the case B, just because of 
the orthogonality of different neutrino mass eigenstates in the flavour space. 
However, they can show up in the case C. There are two ways in which 
the neutrino oscillations might manifest themselves: the first may be called 
global, the second - local. 

To see the global effect one does not need to observe the oscillating term: 
it is sufficient to observe appearance of z/ e , and/or disappearance of v^. This 
could be done without accurate measurements of the time when the decaying 
pion was produced in the target, or the position of its decay point. From the 
theoretical point of view the global effects of appearance and disappearance 
can be described within the approximation of plane waves (monochromatic 
pion). 

To see the local effect, to observe time- and/or space-oscillations, one 
needs an adequate resolution. From the theoretical point of view this also 
requires the pion to be not exactly monochromatic (which is automatically 
the case in any realistic experiment). Otherwise the muons born together 
with V\ would be orthogonal to those born with z/2, and the oscillation term 
would drop out from the total probability. Although our conclusions coincide 
with the usual naive expectations, we feel that it can be useful to present 
this analysis in a little more detail, to avoid any further confusion in the 
literature.^] 

In Part 2 we discuss the decay of a monochromatic pion described by 
a plane wave, derive expression for correlations Pj/^ix^ x v ) (Section 2.1), 
and analyse them in two cases: = (Section 2.2) and 7^ (Section 
2.3). Then, in Section 2.4, we show that if only one of the particles is 
detected, either muon (case B), or neutrino (case C), the oscillating term 

1 Let us note that in the literature neutrino oscillations are treated usually in an 
oversimplified way. The superposition of two neutrino mass eigenstates is described by a 
wave function, not by a density matrix. The two terms of the wave function are usually 
assumed to be plane waves with the same momentum pj, and different energies, E v \ = 
\/pl + ml, E v2 = \/pl + "m\ (see, e.g. books fj, 10, \\ [L2|), or with the same energy E v , 
and different momenta | p>^i | = — m\, | p»^2 1 = \j E% — m\ (see ||, [l3|], and especially 
|l4| , |l5| , |l6|). The 4-momentum conservation in the decay 7r — > \iv is usually ignored. 
Besides U |), there were only a few papers jl7|, |l9), in which possible kinematic 
manifestations of the 4-momentum conservation had been discussed. The neutrino wave 
packet had been considered only by few authors |^o|, ^l], p2| |. 

After completion of this paper we have read section 5.2 in rcf . p"2| , where momentum 
conservation and wave packets were discussed in a spirit close to ours. 
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disappears. We argue that for muons the absence of oscillations is natural, 
but for neutrinos it is an artifact of the plane wave approximation for the pion. 
We also show that the global effects of z/ e -appearance and z/^-disappearance 
are reproduced in the plane wave approximation. 

In Part 3 we use pion wave packet and derive expressions for fiu system, 
as well as for separate beams of muons and neutrinos (Section 3.1). These 
expressions are extremely simple, as the oscillating terms depend only on the 
time of flight of neutrino from its birth till its detection. In Section 3.2 we 
rederive some of the results of Section 3.1 by using the technique of plane 
waves and of the so called "classical a posteriori trajectories" for muons and 
neutrinos. Part 4 contains a brief summary and a few concluding remarks. 

2 Decay of a monochromatic pion state 

2.1 Experiment of the type A: The probability formula 

To begin with, let us assume that pion has a definite 4-momentum p n = 
(E n ,p n ). Then the 4-momenta p^ n = (£ Mn ,p Mn ) and p un = (E un ,p un ) are 
determined by the conservation law, 



and the direction of, say, p un . In what follows we assume that the momenta of 
all three particles are known: either measured, or deduced from kinematics. 
All particles are on the mass shell: 



Ptt Pan ~r Pun 5 



71=1,2, 



(4) 



2222 22 22 2 

Pw = m n, Pul = P u2 = "V Pvl = m l> Pv2 = m 2 ■ 



(5) 



The wave function of the \iv system evolves in space-time as 



+ \!j)e- ip ^- Xi) \u2)e- ip ^- Xi) sm9 = 
= l/^di^e-^ 1 cos#+ \v 2 )e- ilp2 sm6) , 



(6) 



where 




(7) 



<f2 = Pv2{Xv - Xi) + Pn2(Xu ~ X t ) . 
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\v%) and \u 2 ) are "ket" state vectors of the muon p and neutrino mass 
eigenstates v\ and 1/2 respectively; p v \ and p,, 2 are 4-momenta of V\ and z/ 2 , 
and p^i and p M2 are those of muons emitted together with v x and v 2 \ x M 
and x v are space-time coordinates of the muon and neutrino, while is the 
coordinate of the decay point [i for initial). 

The amplitude to detect a muon at a space-time point x M together with 
a neutrino of flavour a (a = e, /i) at the point x u is 

au a {Xf,,x v ) = {yL]P a \ijj p ^{x^x v )) = (3i a e~ lipi + p 2a e~ llf2 , (8) 

and the corresponding probability is: 

■^Va ^^1 Xv ) \ a va{ X Hl X v)\ (^) 
= flla + fila + Zfilafoa COs(y>i - ^ 2 ) , 

where 

<Pi ~ ¥2 = Pvifai/ - &i) - Xi) -prt(x v - ar<) - JV 2 (x M - x<) = (10) 

= (Pki -P^X^ - + (p M i -p^X^/i - = 

= {Pul - Pu2){Xy - Xi - Xf, + Xi) = {p vl - p v2 ){x v - Xp) . 

We made use of the equality p ul — p v2 = — — p^), which follows from 
the 4-momentum conservation (f|). Note that <pi — (p 2 is explicitly Lorentz 
invariant. 

It follows from (|3|), that 

Pip — c °s 2 0> fan = sin 2 9, [3 2e = — flu — sin 9 cos . (11) 
Hence if the neutrino is detected as f M , 

HfK x ^ x ») = cos4 e + sin4 # + 2 sin2 9 cos2 cos (^i - fz) ■ ( 12 ) 
If the neutrino is detected as u e , 

P { v f\x^x v ) = 2 sin 2 cos 2 - 2 sin 2 cos 2 # cos(<^ - <p 2 ) . (13) 

The probability (§) oscillates in space and time with the change of x^ 
and/or x u , presenting a kind of EPR effect. Actually these oscillations de- 
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pend only on the differences Xj, — and t v — t M , and can be observed in an 
experiment, which detects both muon and neutrino from the same decay.0 

The expression (^) was in fact derived but misinterpreted in ref . [|TJ . An 
extra e _r,li in [|lj, which takes into account the decay of muon, is not essential 
for our analysis. 

Before proceeding to the discussion of experiments B and C, let us look 
at eq.(||) a little closer by substituting the expressions for and p u in two 
cases: = and p n ^ 0, the former being a limiting case of the latter 
one. We start with absolutely non-realistic case of vanishing p^, because it 
had been discussed in because it has its own subtleties, and because the 
pion rest-frame values of the muon and neutrino energies enter some of the 
expressions for the more general case. 



2.2 Experiment of the type A: pion with momentum 
strictly equal to zero 

In this situation all four particles (tt, //, v\, v-i) are described by plane waves 
and are fully non-localized. The wave function of the \xv system evolves in 
space-time according to eqs.@ and ([?]). 
In the rest frame of the pion: 

ml-ml + m 2 n 
E ™ = 2^ ' (M) 



ml + m 2 ^ — m 2 



Kn = " , (15) 



2m 



2 In fact, an experiment of this kind has been started at IHEP (Protvino) by S.P.Denisov 
et al. |^3|, |24|. They are looking for (semi)leptonic decays of kaons: K — > fj,v, K — > [ivn 
and K — > evn. The energy of kaons is Ek = 35 ± 1 GeV. The experiment has two 
detectors. The first detector measures muons and electrons coming directly from the kaon 
decay; the second detector measures muons and electrons produced by neutrinos. Thus, 
by detecting a neutrino interaction they are able to determine the momentum of electron 
or muon, which was born together with that neutrino in kaon decay. In this way neutrinos 
are "tagged". By analysing the kinematics they are able to reconstruct position of the 
kaon decay point with an accuracy of 3 m. The tagging provides accurate information not 
only on the position of the origin of neutrino, but also on its original flavour. This allows, 
in principle, a more accurate determination of the neutrino mixing angle 9. 
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P"n = -PL ■ (16) 

After the neutrino and the muon are detected we may conclude that pion 
has decayed somewhere on the line connecting two detection points, but a 
priori we are in principle unable to indicate the position of the pion decay 
point on this line. 

As for the frequency of oscillations and the oscillation length, they are 
determined from eqs.(14)-(|l6|): 



Ap° = P°i - P° 2 , 



v ° Ap °~-^^ = -z^ = A ^-i5' (18) 

where L° = 2E°/(m2 — m|) is the standard oscillation length. Note, that 
AE® and |Ap°| are not equal to each other and to 1/L°. Thus, we have in 
this case in eq.(^) 

E° E° 

^-^m^^Ta^' (19) 

where d is the distance between two detectors (detection points). Taking d 
to be fixed, and by measuring At = t v — t M , one can deduce a posteriori 
the point on the line connecting two detectors, where pion has decayed (see 
Section 3.2). 

2.3 Experiment of the type A: pion in flight with strictly 
fixed (sharp) momentum 

Let us consider now a beam of pions moving from left to right along the line 
which connects the muon and neutrino detectors. f\ 

3 The muon detector has to be to the left of the neutrino detector because the latter 
must be shielded. 
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If the velocity of a pion = p n /E n (v n = \v n \) is low enough, v n < v®, 
where v° = p^/ E®, both /i and v will be detected if pion decays "inside", 
between two detectors. For v n > v® both /i and v will be detected if pion 
decays "outside", to the left of the muon detector. For = u° pion must 
decay in the muon detector. Let us express ipi — ip 2 through the time interval 
At = t u — t,j, and the distance between the detectors d = x„ — x^: 

(P! - v? 2 = AE„At - dAp u = {E vl - E u2 )At - (p ul - p v2 )d . (20) 

Lorentz transformations give 

E un = (E° un + ^p° un h, , P w = (pl + v,0, , (21) 
E m = {E° + v w p° ) 7jr , p„ n = (p° + vAh* , (22) 



where j % = — v\ = E^/m^, and thusf] 

Vl ~ V2 = zS^ ~ «**5) A * + {E ° - v ^ )d ) ■ (23) 

If we take into account that 

E° 

L ° = L t' (24) 

we see that the oscillation frequency is proportional to E-^E^/L, while the 
oscillation length decreases as L/E n E u . The explanation of this drastic de- 
pendence which looks quite unexpected will be given in Section 3.2 



2.4 Experiments of the types B and C with monochro- 
matic pions 

Oscillations described by equations similar to eqs.(P)-(|T3|) and ([19]) were 
referred to in ref.[|IJ as oscillations of a muon beam. In order to clarify 
the situation let us consider first the detection of muons (without detecting 
neutrinos) in the case v n = 0. The beam of muons is described not by a wave 



function, but by a density matrix (see, e.g. ||25[|). The probability of the muon 



4 Let us note that the expressions for AE U and Ap y , as given by eqs.(pO|) and (|2 
differ from eqs.(5.8) and (5.9) for the same observables in ref.flj]. 
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detection is obtained by integrating eq.(^) over the neutrino position Xj, and 
summing over neutrino flavours. Each of these operations results in the 
vanishing of the oscillating term in eq.(Q). Integration over x„ simply leaves 
no dependence on x M . Summing over v a [a = e, //) also eliminates oscillations 
since for the case of final and for the case of final v e the oscillating term 
enters with the opposite signs, see eqs.(|T2"D, ([13]). Thus in the case B the 
probability of detecting a muon does not depend on x M . Moreover, it does 
not depend on the muon momentum p M . In the case = there is no beam 
of muons: they come to the muon detector isotropically because the decaying 
pions are fully de-localized due to uncertainty relation. 

Note that the same refers to the case C: although the global effects of 
z/ e -appearance and z/^-disappearance are present, the flux of neutrinos after 
integration over x M is isotropic, and the oscillating term is washed out. The 
latter conclusion may look to be in contradiction with the generally accepted 
theory of neutrino oscillations. However the problem is solved as soon as we 
take into account that in real experiments we never have a decaying pion with 
sharp momentum. Its momentum distribution has a non-vanishing coherent 
spread. As a result pions are described by coherent packets of plane waves, 
and are localized in space-time, [25fl- p9[l . 

Let us now repeat this reasoning in the more formal way; namely, come 
back to the probability formula (|9|) and use it for formal description of ex- 
periments of the types B and C in the general case of a plane wave with 
sharp momentum p^. This is straightforward: one should just sum over all 
the states of neutrino and muon, respectively. 

We begin with the observation of muons (case B). The states of neutrino 
are labeled by the "flavour" index a = e,fi and by the "position" Xj,. Thus 
the probability of the \x detection at the space-time point x^, is equal to 

P^M = E / P£X*» = (25) 



ebc. 



This can be alternatively formulated as a transition from the wave function 
of \iv system to the density matrix ]25[ for \i given by 



Pyix^x'^) = Tr u / i[)(x fM ,x 1/ )'4)*(x' fl ,x v )d-x u . (26) 
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Eq.(|25|) is identical to fl2"6|) for x M = x' . In accordance with the general 
principles of quantum mechanics, the r.h.s. of (25) and (^) are automatically 
independent of t u . 

It remains to substitute (|9[) into (|25|). Integration over gives rise to 
8 (p Ul — p U2 ) = 0, because ^ p U2 . This is already enough to eliminate 
the oscillation term. However, it vanishes in experiment of the type B for a 
more fundamental reason: because of the summation over neutrino species. 
Indeed, it is easy to see that the oscillation term cancels in the sum of (|12]) 
and (|l^) even before integration over x^. Thus we get 

Pi B) M = 1 • (27) 

An attentive reader can get suspicious at this point. We identified two 
reasons for elimination of oscillations in the case B: orthogonality of neutrino 
the states in the flavour and momentum spaces. If we now turn to the 
observation of neutrinos (case C), the first reason is absent (a muon is always 
the same particle), but the second reason is still present: p^i ^ p At2 - and this 
is enough to eliminate the neutrino oscillating term in the case C, contrary 
to any reasonable expectations. 

Still, we insist that in the "sharp" case under consideration this is really 
true: the probability to observe neutrino v a at point x v 

Pl c a \x u )= Jp^\x^x u )d^ (28) 

does not contain oscillating term if one substitutes (^) into (|28|). If one uses 
eq. ([12]) for P^ A '{x^ x v ) and eq. ([l3|) for P^i^x^, x u ), one gets: 

PiTM = Pi + = cos4 9 + sin4 9 < 1 > ( 29 ) 



P<& ( Xv ) = (3f e + (3l e = 2 sin 2 6 cos 2 9 > . (30) 

Thus, global manifestations of the oscillations, the appearance of u e and 
disappearance of z/ M , are evident, but the oscillating terms themselves are 
absent. 

There is a simple reason for this: our assumption that the decaying pion 
has a definite momentum, which is never true in experiments. It is more 
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or less obvious that allowing a minor coherent dispersion in the momentum 
distribution of pion, one will always find solutions to the equation 



pAp*) = pAp*) • (31) 

As we will see (footnote 8), eq.(|31|) has a solution at l^p^l = \p n — p n \ ~ 
1/(1 -v v )L. 

The momentum dispersion is necessarily present in any realistic exper- 
iment: pions are usually localized in a region much smaller than l/|5p„.|, 
and then, by the uncertainty relation, their momentum dispersion^ should 
be much larger than j^p^. 

Thus we need to return to the very beginning and repeat our analysis, 
allowing some small (as compared to the masses and energies in the prob- 
lem), but non- vanishing momentum dispersion in the wave function of the 
original pion. This is a simple calculation, but still it deserves being done: 
for example, one should check that such dispersion does not wash out the 
oscillations in experiments of the types A and C. We shall also use this new 
calculation to represent eq.(||) in a somewhat different form; this can help 
one to better understand the "paradoxical" results of Section [275. 



3 Decay of a pion with small momentum dis- 
persion 

3.1 Experiments A, B, and C 

Let us now assume that pion has been created at a space-time point with 
some momentum distribution <f)(p n ). This means that at the decay point Xi 
the pion wave function is 

~ O = / 0(p.)e-^-^p., vl = ml, (32) 

To avoid confusion, let us emphasize that here we speak about the coherent dispersion 
of a pion produced in a given act of collision with accompanying particles being in a given 
state. This should not be mixed up with non-coherent momentum dispersion of pions in the 
same beam. Non-coherent means that the pion is produced with the same accompanying 
particles, but being in a different state, or with other accompanying particles, or produced 
in a different collision act. 
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and that of the emerging \w system, 

V{x ll ,x v ) = J <j>{p v )e^( Xi -^ ) 4>j ht {x ll ,x 1 ,\x i )dp ir , (33) 

with ip Pl ,{x^x v \xi) given by eq.(|). 

Actually ^(x^, x„) does not depend on x%, because we imposed the con- 
servation law (f|) in all our formulas. 6 We shall, however, proceed a little 
differently and keep for a while. Moreover, we introduce the condensed 

By using eq.([33|) we define the amplitude A Ua (x fl , x v ) (compare with eq.(§) 
for a Va {x^x v ))\ 

^fal-^A 4 ; •Ev) = I < ^P7r0(P'7r) ^ ] ftna 6Xp ( Vp^Xfa iPfMn^fii ^Vvrv^vi) > (34) 

and get for the probability: 

^ ) (^,^) = |A a (x M ,x ! ,)| 2 = (35) 
= y /" rfp 7r rfp 7r 0(p 7r )0(p 7r )^/3 na /3 fla e" 



-ii'fin—fn) 



where 

and n(n) = 1, 2. 

Let us now make use of the assumption that dispersion of the distribution 
4>{p n ) is small as compared to the typical energies in the problem. This allows 
one to put 

5E n =E n -Eir= -^tPtt(Ptt ~ Ptt) = V 7r (p 7r - p^) = V n 5p n , 
5E^ = E^ n - E^ n = — P^(p M „ - p^fi) = V /t (p M n - P^n) = V^p^ , (37) 

Ef, 



6 The usual logic in quantum field theory is reversed: one includes integration J d 4 Xi 
in the definition ( p3| ) of ^{x^^Xv), and this integration leads to the conservation law 
Our presentation is organized to minimize the length of the reasoning. We do not need 
to be too careful about such details; note that we ignore all the normalization factors, 
irrelevant for the main topic of the paper. 
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2 2 

mi — mi 



5E V = E un — E un = — \ T> u {$i vn — p u fi) H ^ — - 



2 2 2 2 

mi — mi „ m„ — 



v„(p™ - p^s) H ttt; — - = v„5p„ + 



2E„ r " 2E tJ 



where v's are velocities of the particles, v = p/E. (Formulas fl3~7|) are derived 
by subtraction of the mass-shell equalities E 2 = p 2 + m 2 and E 2 = p 2 + m 2 ). 
Note also that from energy conservation one gets 

5E n = 5E n + 5E U . (38) 

Substituting this into fl36|) we get: 

<fn - <fn = — (Pit _ P7r)(x i7r - V„-£ ilr ) - 

-(P/m - P M n)(x M j - v^) - (p^ n - p,J(x^ - v v t vi ) + (39) 



2E, 



V 



From eq.(^) we may derive formulas for the cases A,B,C with the pion 
being described by a wave packet. 

When considering the case A, let us note that the first three terms in 
eq.(|3"§|) vanish on the trajectories of the particles, 



H7T 



HenceQ 



- y?n = -j- for n = 1, n = 2 . (41) 



Note that in the case A the muon detector is used in order to deduce the 
position in space-time of the pion decay point i. Hence the value of t v i (in 
the case A we may call it tA) can be determined without the knowledge of 



the moment of the pion production t w (x n = (i^x^), see eq.(|32|)). This will 



7 Let us check again that the phase difference given by eqs.(^)-(^T]) is Lorentz invariant. 
In the reference frame, which moves with the velocity u, t V i — > ~f u {tvi + ux„.;) = (1 + 
v^u)lutui and E v — ► J U {E V + upj = (1 + VlV v )^ u E v , where j u = - u 2 , so that 

t U i/L — {tyi/ E v ){m\ — m\)/2 remains invariant. 
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be done explicitly in eqs.(|55|)-(|58|). Our new expression for the probability 
in the case A, 



2 1 

mi — mi 



' Xu ) ftnaftna 6 X P ( i ^j^, t p 



= Pla + PL + Wlafoa COS ^ (42) 

looks absolutely different from the old one (compare (^) and ([42]) with (|9|), 
(TU), (|T^), (^))- In the next Section we will show that in fact they are 
equivalent, and ipi — ip 2 = <fi — y?2 = t ui /L. 

The same momentum dispersion technique can be applied to the descrip- 
tion of realistic experiments of the types B and C . In the case C (sin- 
gle neutrino detector) by integrating (|35|) over x M , one gets delta- function 
$(Pfm ~ P/xn)) where p Mn and are expressed through p^ and p„-, respec- 
tively If n ^ n, the argument of the delta-function vanishes for non- vanishing 
Ptt — Ptt = 8pn, which for high energy pions must be equal to 2^1 / L = 'Jx/L 
(see also eq.(|23|))f| The corresponding probability is 



Vi c a \x u ) =0Z a + Pl l + 2£l3 la p 2a CM^, (43) 

where 

p = I dpMp^ - 5p w /2)(j)(p 7r + 6p w /2) 

4 /^0 2 (P.) ' 1 } 

and can be expressed through x v and x n , the space-time point where pion 
has been created; in the collinear case (p^Hp^): 

& = d — • ( 45 ) 



8 To find exact expression for p w — p,r = 8pn, let us consider p^ as a function of p w and 
m v , and calculate the difference 5p^ = Sp^ — 5p v . By using eqs.(|37|) and ( |3§| ) it is easy 
to find 5p^, corresponding to 6p M = p M „ — p pR = (and hence SE^ = 0). The condition 
<)Pn — means 5p v = Sp^ and 8E V — SE^, hence (v ff — v^p^ = 1/L, where we put 
5ml — m\ — m\. For the collinear case, v„ || v f , we get |Jp ff | = L _1 /(l — u^) since i>„ w 1. 
For high energy pions |5p w | is much larger than L~ 1m . \6p„\ ~ = ^/L . 
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For monochromatic pions £ = 0, and 



Pj?(**l£ = 0)=/& + & ) w 
while in realistic experiments £ = 1, and 

P£> (*„ | £ = 1) = & + dla + Wlafoa COS f . (47) 



The value of in eq.(fl^) can be fixed only if the moment of the pro- 
duction of a pion in target, t n , is known with a very high accuracy. For that 
purpose the time-structure of the proton beam should be of the order of few 
nanoseconds. Otherwise the oscillating term in eq.(f43"D would be washed out, 
and only global manifestations of neutrino oscillations ( ^-appearance and 
z/^-disappearance) would remain; compare with eqs.(p9|) and (|30|). Another 
obvious way to fix and thus to observe the oscillating term is to ascertain 
the position of the pion decay point Xj. 

As for the case B, it has been already mentioned that the oscillating 
term vanishes after summation over neutrino flavours (see Section 2.4). This 
statement is obviously true in the situation when pion is described by a packet 
of plane waves: summation over a = e, /i using eq. 

(|42|) gives V^ B \x^) = 1. 



3.2 On equivalence of the two representations for P^ A > 



There are two main ingredients in the previous section, namely: 1) the mo- 
mentum representation of the pion wave packet (eq.([!3|)), and 2) the classical 
trajectories for fast particles, eq . (|40|) . Actually we could use only the second 
one in order to get eq.(fE^), and to solve the "paradoxes" of sections 2.2 and 
2.3. 

In this section we will work with plane waves and at the same time with 
classical trajectories for /i and v (see e.g. f26|l ). After \i and v are detected 
at certain space-time "points" with given momenta, we have, so to say, "a 
'posteriori packets" of these particles, for which x Mi = v^t^, x. ui = v u t ui , 
where v M = p^/E^, and v„ = p u /E u . Note that the wave length of \i or 
v does not exceed 10~ 13 cm, while the characteristic distances in neutrino 
experiments are larger than hundred meters. 
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For v and /i on the mass-shell: 



2 2 

m\ — m\ 



AE V = E vl - E u2 = v„Ap v + 1 1 , (48) 

l tii,, 



J v 



AE^ = - E^ = v^Ap^ , (49) 

where 

Ap v = p ul - p u2 , Ap^ = p Ml - p m2 • (50) 
Note that if the pion momentum is sharp, 

AE W = , Ap^ = , (51) 

then 

AE^ = -AE V , Ap M = -Ap„ , (52) 

and hence 

v.Ap, + v M Ap M = (v„ - v„)Ap„ = — ^ = - . (53) 

In eq.fl4"T|) we have shown that tp\ — (p^ = t ui /L. Consider now (pi — (p 2 as 
given by eq.(|lOD, and use eqs. fl48|) , (fSf), and (|52f) : 

</>l - V?2 = (Pul - Vu2)(Xvi - X,j,i) = 

= AE u {t vi - t^) - Ap F (x„ - x M i) = (54) 
= t H (AE v - v u A Pu ) + t^AE^ - v,A P/1 ) = j- . 

Thus <fi — (p~2 = V^i — ^2) as was promised after eq. ([42]) . Eq.(|54|) directly 
brings us from eq.(|D to eq . (|42|) . We would get the same result in the cases 
of a pion with p^ = (eq. flT^)), and of a relativistic pion (eqs.(|2"0"|)-(P3"|)). 
At fixed distance between the detectors d the measurement of At (time 
difference between "clicks" of the two detectors) allows one to find the space- 
time point Xi of the pion decay. Thus for a pion with v„ > t>° decaying to the 
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left of the muon detector we have in the collinear case discussed in Section 
2.3 



\^-(t Xj 



v^-ti), |xj,-Xj| =v v (t v -ti) , (55) 
hence 



(56) 



and 



tt t = t vi = t v -U = v ^ ~ U) + K ~ ^ = ~ V " At + ■ (57) 



Taking into account that v v ~ 1, we see that in the last equation denom- 
inators are extremely small for ultrarelativistic muons. That means that the 
decay of the pion discussed in Section 2.3 must take place at very large dis- 
tance from the detectors, unless the numerator in eq. (|57D is also very small. 

For a pion with v n < v®, decaying between the two detectors 



, , v u At + d 



~ Ui ~ ^+77 ' (58) 

It is important to stress that in eq. (p9|) or fl54[) the phase <p± — <y2 2 depends 



only on t vi and does not depend on t^. Thus these formulas simply describe 
the standard neutrino oscillations from the point of creation of the neutrino 
till the point of its detection. Accordingly they do not depend on the position 
of the muon detector. Therefore EPR-like correlations between muon and 
neutrino detection appear only when we express t ui in terms of At = t u — t M 
and d = — x M . If we assume that a special detector measures the decay 
point of the pion Xi, then the situation becomes absolutely trivial. 



4 Concluding remarks 

Let us summarize the main results. 

In consideration of decay ir — > fiu and subsequent neutrino oscillations 
we assumed that the 4-momentum is conserved and all particles (n, /i, z/'s) 
are on the mass-shell. 



17 



We used two different approximations: 

1. The pion is monochromatic (has definite, sharp momentum). Then the 
same is true for \i and z/'s (up to an ambiguity in rotation of the \iv plane, 
which is inessential for our purposes). Thus all particles are described by 
plane waves. Non-trivial phenomena, associated with neutrino oscillations, 
are due to a slight difference between the 4-momenta of neutrinos with dif- 
ferent masses. As a corollary, the muon, created in decay of the same pion 
together with different neutrinos, also possesses slightly different 4-momenta. 

2. The pion is localized in space and is described by a (coherent) wave 
packet. Then the \w wave function is a linear superposition of wave packets, 
resulting from decays of pion at different times. It is reduced to a product of 
packets for ji and u, when any of the particles is detected at a definite space- 
time point - this gives rise to what we call "classical a posteriori trajectories" . 

We analyzed the following three experimental situations. 

A. Muon and neutrino detected in coincidence. The probability of de- 
tecting fM at a space-time point x^ and u a (a = e, ji) at a point x v is given by 
eq. (|4"2|) , where is defined by eq.(|57|) or (|58|). 

B. Neutrino is ignored, while fi is detected at the space-time point x M . 
The corresponding probability does not depend on x M at all, V^ B '(x,j) = 1, 
and no traces of neutrino oscillations are seen. 

C . Muon is ignored, while v a is detected at the space-time point x v . This 
time the probability depends essentially on the spread of the (unobserved) 
muon wave packet and thus on that of original pion. For monochromatic pion 
or, more precisely, for |5p„-| C 7^/1 in a collinear case, the parameter £ in 
eq. (|43|) vanishes, and there is no oscillating x v dependence in the probability, 
see eq.([|6]). The only global consequence of neutrino oscillation in this case is 
that V^p > 0, while Vjf^ < 1. For large enough[] dispersion, l^p^l ^> 7^/£, 
we have £ = 1, and the standard expression V^'{x u ) given by eq.([I7|). If 
tjr is not specified (not known in a given experiment), the probability should 
be averaged over t n . This eliminates the oscillation term, but preserves the 
global dependence on neutrino flavour v a . Note, however, that t n was used 
by us to deduce the position of the pion decay point x«. If this position is 
known from other considerations (e.g. short decay pipe), then t vi is fixed, 
and the oscillating term is of course retained. 

9 Note that for ultrarelativistic pions, and in the collinear approximation, ^p^ has to 
be many orders of magnitude larger than the naive estimate, l^p^l S> 1/L ! 
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When this paper was near its final stage of preparation, the authors 
of ref. KJ have replaced the original paper (E.Sassaroli, Y.N.Srivastava and 
A.Widom "Charged Lepton Oscillations" , |hep-ph/9509261| ; September 1995, 
15 pages) by another one: Y.N.Srivastava, A.Widom and E.Sassaroli "Lep- 
ton Oscillations", |hep-ph/950926T| v2, (November 24, 1996; 2 pages) with the 
same figures and the same claim formulated in a new abstract: "A simple but 
general proof is presented to show that Lorenz covariance and 4-momentum 
conservation alone are sufficient to obtain muon oscillations in pion decay if 
the recoiling neutrino oscillate." 

On December 13 a new paper |30|| by A.Widom and Y.N.Srivast ava has 
appeared in the hep-ph Archive. In this paper they claim that muon oscil- 
lations associated with mixed neutrino mass matrices should manifest them- 
selves in the experiment measuring (g—2) - the anomalous magnetic moment 
of the muon. As is known, the (g — 2) experiments are not done in coinci- 
dence with detection of neutrinos which accompany production of the muons. 
Therefore we believe that no effect associated with neutrino mixing should 
be seen in (g — 2) experiments. 

As follows from our discussion, oscillations of correlation probability may 
be observed in the two-detector experiments which measure both charged 
leptons and neutrinos in coincidence. Q These oscillations would look like 
an EPR effect, though this is a simple consequence of standard neutrino 
oscillations and relativistic kinematic relations. Muons do not oscillate. 



10 One is tempted to apply arguments developed in this paper not only to hypo- 
thetical neutrino oscillations, but also to the already observed neutral kaon oscillations 
and to similar oscillations of neutral heavy mesons D°, D Q S1 B Q and B° s . The obvi- 
ous analogs of the it — * [iv decay would be the decays A+ — > pK°, S° — > AK°, or 
D+ -> K+K , D+ -► if°7r+, D+ -> K+K°, or decays of B-mesons: B+ -> D°ir+, B+ -► 
D Q D + , B + — ► D°D+ . The corresponding decay widths in the above examples are larger 
than the oscillation frequencies, and therefore the dispersion of momenta is provided by 
these decay widths. The description of EPR effect in <t> and T decays in terms of amplitudes 



has been advocated in ref. 31 



Claims that neutral kaon oscillations must produce oscillations of A-hyperons, and that 
the frequency of kaon oscillations in 0-meson decay is factor of 2 larger than its standard 
generally accepted value, have been published in refs.|32| [53|, and have been disproved in 
refs.fl3l M. 
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